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Abstract 



^ We derive a conservative coincidence time window for joint searches of gravita- 

2 tional-wave (GW) transients and high-energy neutrinos (HENs, with energies 

^ > lOOGeV), emitted by gamma-ray bursts (GRBs). The last are among the 

>• most interesting astrophysical sources for coincident detections with current 

^ and near-future detectors. We take into account a broad range of emission 

^ mechanisms. We take the upper limit of GRB durations as the 95% quantile 

of the Tgo's of GRBs observed by BATSE, obtaining a GRB duration upper 
limit of ~ 150s. Using published results on high-energy (> lOOMeV) photon 
light curves for 8 GRBs detected by Fermi LAT, we verify that most high- 
energy photons are expected to be observed within the first ~ 150s of the 
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GRB. Taking into account the breakout-time of the relativistic jet produced 
by the central engine, we allow GW and HEN emission to begin up to 100s 
before the onset of observable gamma photon production. Using published 
precursor time differences, we calculate a time upper bound for precursor 
activity, obtaining that 95% of precursors occur within ~ 250s prior to the 
onset of the GRB. Taking the above different processes into account, we ar- 
rive at a time window of tnEN — tow ^ [— 500s, +500s]. Considering the 
above processes, an upper bound can also be determined for the expected 
time window of GW and/or HEN signals coincident with a detected GRB, 

tew - tcRB ~ tHEN - tcRB & [-350s, -M50s]. 

Keywords: 
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1. Introduction 

The joint detection of gravitational waves (GWs) and high-energy neutri- 
nos (HENs) from astrophysical sources in a multimessenger approach promises 
significant scientific benefit over single messenger searches. Such advantages 
include increased sensitivity, as well as additional information on the un- 
derlying mechanisms responsible for creating such multimessenger sources 
[1, 2, 3, 4, 5]. 

GW and HEN detectors have made significant progress in recent years. 
The first generation of GW observatories, such as LIGO [6], Virgo [7] and 
GEO [8] have already made astrophysically interesting observations (e.g. [9, 
10, 11, 12, 13, 14, 15]). The second generation of these detectors (e.g. [16]) 
is expected to be roughly an order of magnitude more sensitive around ~ 
150Hz. HEN observatories currently under operation are IceCube [17], a 
cubic kilometer detector at the geographic South Pole, and ANTARES [18] 
in the Mejditerranean sea. IceCube has reached its design volume with 86 
strings, while ANTARES is operating at its design volume with 12 strings. 
ANTARES is planned to be upgraded into a cubic kilometer detector called 
KM3NET in the following years [19]. The Lake Baikal neutrino detector is 
also planned to be upgraded to a cubic kilometer volume aiming to detect 
HENs [20]. 

One of the most promising class of sources for joint GW+HEN searches is 
represented by gamma-ray bursts (GRBs). GRBs are thought to be produced 
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by internal shocks in relativistic expanding fireballs from energetic cosmic 
explosions [21]. Fireballs [22], i.e. the relativistic outflows or jets of plasma, 
are created when the central engine releases a large amount of energy over a 
short time and small volume [23]. These jets are often thought to be launched 
along the rotational axis of the progenitor, powered by the gravitational 
energy released during temporary mass accretion onto the central black hole 
[24]. 

There are two main cosmic progenitor candidates that are expected to 
produce the relativistic outflows and shocks of magnetized plasma described 
above. Long-duration GRBs (>2s) are thought to be produced by stellar 
core-collapse [25], while short GRBs (<2s) are predominantly associated with 
compact stellar mergers [26]. 

Joint GW+HEN searches are of special importance for sources dark in 
other messengers, e.g. that have little or no electromagnetic emission. Such 
potential sources include choked GRBs [27]. If a stellar core-collapse emits a 
relativistic jet from its core that emits gamma-rays, these gamma-rays will 
only be observable from the outside once the relativistic jet broke out of 
the stellar envelope. Relativistic jets that stall before reaching the surface 
of the star can still emit observable HENs [27, 28], while the envelope is 
transparent to GWs. The result is GW+HEN emission with little or no 
gamma-ray counterpart. 

Previous single or multimessenger searches for GWs, HENs and/or elec- 
tromagnetically detected GRBs have defined various different time windows. 
Aso et al. [2] designed a GW+HEN multimessenger search algorithm and 
investigated possible time windows from 0.1s to up to 1 day. In a search 
for HEN counterparts of detected GRBs, Abbasi et al. (IceCube Collabo- 
ration) [29] consider three different time windows around each burst. The 
first time window covers the observed prompt gamma-ray emission of the 
GRB, expecting prompt neutrino emission to overlap with prompt gamma- 
ray emission. The second considered time window covers neutrino-producing 
processes before gamma-rays can escape from the fireball. This time win- 
dow is taken to be 100s immediately preceding the prompt emission time 
window. Abbasi et al. (IceCube Collaboration) [29] also consider a generic 
time window that is chosen to be much wider than the previous two to in- 
clude possible unknown mechanisms. This third time window is chosen to 
be [— Ih, 3h] around the start of the prompt emission, a compromise between 
including unknown mechanisms and keeping the HEN background low. 

In a HEN search for GRB 0803 19B, one of the brightest GRBs ever ob- 
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served, the IceCube Collaboration [30] uses two time windows for HEN de- 
tection. A shorter time window of 66s was used that overlapped with the 
GRB's prompt gamma-ray emission. Another, ~ 300s long extended time 
window was also analyzed. Data from the detector was only available in this 
300s around the GRB which motivated the choice of this time window. 

Very recently a search for neutrinos from 117 GRBs was conducted using 
data from the IceCube detector's 40 string configuration [31]. Two different 
searches were performed; one model independent search, using no model 
for the energy distribution of HENs from GRBs. This search used a time 
window of ±24 hours around the burst. The second, model dependent search 
used the predicted energy distribution from GRBs by Guetta et al. [32]. 
This search considered the observed start {Tstart) and stop {Tstop) times of 
gamma-ray emission from each GRB. A probability distribution was assigned 
for expected neutrino times with uniform probability between Tstart and Tgtap, 
and with Gaussian tails around this uniform window. The width of the tails 
was chosen to be Tgtop — Tstart, and was constrained to at minimum 2 seconds 
and at maximum 30 [31]. 

In GW burst searches triggered by long and/or short GRBs, performed 
with the LIGO- Virgo GW network, [33, 34, 13, 35], a time window of [-120s, +60s] 
was used around the GRB trigger. This time window is longer than the dura- 
tion of most analyzed GRBs. It was also chosen to include models predicting 
gamma-ray emission up to 100s after initial GW emission (e.g. [36]), as well 
as measurement uncertainties in GRB trigger time. For a search for GW 
inspiral counterparts of detected short GRBs, Abadie et al. (LIGO Scientific 
Collaboration and Virgo collaboration) [37] used a time window of [—5s, -fls) 
around the trigger time of GRBs, aiming to capture "the physical model with 
some tolerance for its uncertainties." 

The present article, for the first time, discusses the temporal structure^ 
of GW and HEN emission from GRBs relative to each other, with the main 
goal of defining a common on-source time window. As many aspects of GRB 
emission processes are still debated, we base our analysis largely on model- 
motivated comparisons with observations. 

The article is organized as follows. We first (Section 2) present a summary 
of the obtained GW+HEN coincidence time window building on the results 
and discussion of Section 3. In Section 3 we review the potential emission 
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processes in detail, using model-motivated comparisons with observations 
and simulations to constrain the duration of each emission epoch. We refer 
the reader to Section 2 for an overview and summary of the results, while 
the details can be found in Section 3. 

2. Summary - Time Window of GW+HEN Emission 

In this section we determine an upper bound for the time difference of 
HEN and GW signals from GRBs. We will refer to this time interval as 
GW+HEN coincidence time window. We aim to obtain an upper bound 
in which HENs and GW transients are most likely to arrive. The section 
summarizes the expected durations of different emission processes that are 
discussed in details in the rest of the article. 

The analysis mainly relies on model-motivated comparisons with observa- 
tions and simulations. We consider different types of information we have on 
GRBs from their detected electromagnetic radiation, and infer to the plau- 
sible time frame of GW or HEN emission based on this information. We use 
observations by BATSE [38], Swift [39] as well as Fermi LAT [40]. We esti- 
mate the duration of different emission epochs, and combine them to obtain 
an overall time window. While a joint search for GW and HEN signals mainly 
aims to detect astrophysical sources with unobserved electromagnetic (EM) 
counterpart, we constrain their emission based on the emission of detected 
sources, assuming that these two types have similar emission. 

A summary of the considered emission processes and their time frames 
are shown in Figure 1. Note, that for the case of GW emission, a time 
frame possibly having GW signals means that there can be one or more 
short transients at any time of the time frame. There can also be a longer 
GW emission during this time (see Section 3.3 for further details). 

• We start with the onset of gamma-ray emission from the main GRB. 
The duration of the main GRB is, in 95% of cases, shorter or equal than 
150s (see Section 3.1). During this time frame, as the central engine 
needs to be active, both HEN and GW emission is possible. 

• Recent Fermi/LAT observations (e.g. [41, 42]) have shown that at least 
some GRBs also emit high-energy (GeV) gamma-rays, with temporal 
structure different from that of gamma-rays of typical energies (see 
Section 3.2). Our estimated time scale of the emission shows that most 
(^ 90%) of the high-energy photons are expected to be detected within 
150s after the onset of the GRB. 
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(b) (d)gamma > 100 MeVi 



(e) 500s 



Figure 1: Summary of upper bound of GRB emission process durations taken into account 
in the total GW+HEN coincidence time window, (a) active central engine before the 
relativistic jet has broken out of the stellar envelope; (b) active central engine with the 
relativistic jet broken out of the envelope; (c) delay between the onset of the precursor 
and the main burst; (d) duration corresponding to 90% of GeV photon emission; (e) time 
span of central engine activity. Overall, the considered processes allow for a maximum 
of 500s between the observation of a HEN and a GW transient, setting the coincidence 
time window to [—500s, 500s]. The time window for GW or HEN signals from the onset 
of the GRB is tew — toRB ~ tHEN — tcRB & [—350s, +150s]. Note that we show a period 
between the end of the precursor emission and the start of the main GRB with no GW or 
HEN emission. While we cannot exclude the possibility of GW or HEN emission in this 
period, such emission would have no effect on our estimated time window. The top of the 
figure shows a schematic drawing of a plausible emission scenario. 
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• GW and HEN emission can commence up to 100s before the onset of 
gamma-ray emission (see Section 3.4) due to the time it takes for the 
relativistic jet to break out of the star (e.g, [27]). The central engine is 
active during this time, for the jet needs to be actively driven in order 
to advance. 

• Precursor emission can begin up to fP^^'^'^'^^^"^ ~ 250s prior to the onset 
of the main GRB (see Section 3.5). Since a precursor can be active 
practically any time within this time frame, and since the activity of 
the central engine is a prerequisite of precursor emission, both HEN 
and GW emission is possible during this period. 

• The activity of the central engine responsible for a precursor can begin 
up to 100s before the onset of gamma/X-ray emission by the precursor 
(see Section 3.5). During this time, as the central engine is active, both 
HEN and GW emission is possible. 

• We note that our choice for the GW/HEN coincidence window is not 
affected by what happens in between the precursor(s) and the main 
GRB. Since HEN and/or GWs may be present both at the precursor 
and during the main GRB, the GW/HEN coincidence window needs 
to extend anyway from the precursor up to the end of the main event. 
Thus, such a choice also covers the case in which GWs and/or HEN 
are also emitted between the precursor and the main GRB. 

Considering the processes described above, we obtain a time frame of 500s 
during which both HEN and GW emissions are possible. This means that 
the difference between the arrival time t^EN of a HEN and the arrival time 
tow of a GW transient from the same GRB is within the upper bound (i.e. 
time window) of 

tHEN - tew e [-500s, +500s] (1) 

Note that multiple HEN (or GW) signals from a given source should arrive 
within a 500s time window. 

We emphasize that the above time window is an estimated upper bound 
for GRBs. While the emission mechanisms considered are connected to the 
progenitor associated with long GRBs, short GRBs likely have shorter emis- 
sion span and are therefore also bound by this time window. Also, we exclude 
potential GW/HEN emission from GRB afterglows from the time window 
(the emission of ultra-high energy 10^ — lO^GeV) neutrinos may be more 
extended; see Section 3.4). Other effects, such as some predicted by different 
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Quantum Gravity models (e.g. [43]), can introduce additional time delay 
between different messengers. Searching for such effects may require larger 
time windows [43]. 

The processes considered above can also be used to set an upper bound 
on time window of GW or HEN signals coincident with an observed GRB. 
If the onset of a GRB is detected at tcRB, both GW and HEN signals are 
expected to arrive within the time window 

tew - tcRB = tHEN - tcRB € [-350s, +150s] (2) 

3. Emission Processes in Gamma-ray Bursts 

The widely accepted model of long GRBs, the so-called collapsar model, 
associates long GRBs with the collapse of the core of a massive star, that 
drives an outgoing relativistic outflow (see, e.g, [27]). According to the fire- 
ball model of GRBs [22], gamma-rays, as well as HENs are produced by 
shocks in this relativistic outflow. 

The plasma outflow from the central engine reaches relativistic velocities 
only after leaving the Helium core of the star [27] . Consequently, the time it 
takes for the outflow to break out of the star is approximately equal to the 
duration of crossing the Helium core (tne ~ 10 — 100s) [27] (see also [44]). 
The outflow can advance only as long as it is driven by the central engine, 
therefore the duration of the outflow has to be longer than tne in order for 
the outflow to break out of the star. Due to this requirement, there can be 
stellar core-collapses where the outflow is unable to break through the stellar 
envelope [45]. In these so-called choked GRBs, whose number may exceed the 
number of successful GRBs [45], GWs, low-energy neutrinos and potentially 
HENs can escape the star, while the star remains dark in gamma-rays. 

Besides long GRBs, another distinct subclass of gamma-ray bursts are 
short-duration, hard-spectrum GRBs, with duration < 2s (e.g, [26]). At 
least some of these short GRBs are associated with neutrons star - neuron 
star or neutron star - black hole binary coalescenses (e.g, [26]). Such binary 
coalescence constitutes the central engine that create relativistic outflow as 
discussed above for long GRBs, responsible for the emission of gamma-rays 
and high energy neutrinos. 

This section discusses the different emission processes with special em- 
phasis on their relative timing, defining their contribution to the GW+HEN 
coincidence time window. 
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Figure 2: (left) Distribution of Tgo's for GRBs measured by BATSE (straight line, 1234 
GRBs, [38]) and Swift (dotted line, 534 GRBs, [39]). The vertical line shows tg'^^ , i.e. the 
95% quantile of the Tgo's of BATSE GRBs. (right) Distribution of time difference between 
the onset of a precursor and the onset of the main GRB. The vertical line shows fP'^'^'^'^^''"^ ^ 
i.e. the 95% quantile of the time differences (the delays were obtained by Burlon et al. 
[46] using data from BATSE [38]). 



3.1. Gamma-ray Emission 

The observed temporal structure of gamma-ray emission from GRBs can 
be used as the foundation in the description of the temporal structure of GW 
and HEN emission. For the purposes of GW+HEN analysis, we define a prac- 
tical upper limit for GRB duration as the 95% quantile of the Tgo's of GRBs 
detected by the Burst Alert and Transient Source Experiment (BATSE). We 
obtain this limit using data from the 4th BATSE GRB catalog (1234 GRBs, 
[38]), that includes 1234 GRBs with duration information (see Figure 2). As 
the measure of GRB durations, we use their T90, the time interval in which 
the integrated photon counts from a GRB in the BATSE detector increased 
from 5% to 95%. See Figure 2 for the distribution of Tgo's of BATSE GRBs 
in comparison with tg^^. The obtained duration upper limit is 

tg^^ ^ 150 s (3) 

For comparison we note here that the 90% and 99% quantiles are tffj^^ ~ 100s 
and tg^^ ~ 300s. Considering a similar analysis for GRBs detected by 
the Swift satellite (534 GRBs, [39]), we find that 87% of Swift GRBs have 
T90 < 150s. For Swift, the 90%, 95% and 99% quantiles are 180s, 300s and 
500s, respectively. See Figure 2 for the distribution of Swift Tgo's. 
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3.2. High-energy Gamma-ray Emission 

While typical gamma-ray energies from GRBs are in the range of 0.1 — 
1 MeV, some GRBs also emit gamma-rays with energies above 100 MeV 
(e.g. [41, 42]). Such high-energy component was first observed by the Ener- 
getic Gamma-Ray Experiment Telescope (EGRET, e.g, [47, 48, 49]), and has 
recently been identified for more than 20 GRBs [50] by the Large Area Tele- 
scope (LAT) on the Fermi satellite [40]. GeV gamma-ray emission from these 
GRBs, in most cases, shows radically different characteristics [41, 42, 51] 
compared to MeV gamma-ray emission. 

The different temporal and spectral behavior of GeV photons imply that 
they are likely produced by different processes during a GRB than MeV pho- 
tons. There are different proposed mechanisms responsible for the creation 
of GeV photons. It is also possible that different mechanisms are responsible 
for the GeV photon output for different GRBs. Studying the high energy 
emission of GRB 090510, Ghirlanda et al. [52] suggested that this may be 
explained entirely due to external shock emission. On the other hand, several 
other authors [53, 54, 55] have shown that internal shocks may also play a 
role, especially during the first few seconds of the high-energy tail. Toma et 
al. [56] argue that GeV photon emission can be due to Compton upscattered 
photospheric emission, during which photons are upscattered by electrons in 
internal shocks. Photohadronic interactions in proton-dominated bursts have 
also been considered, that may be responsible for the production of the high- 
energy component [57, 58]. In this case, GRBs with GeV photon emission 
can be excellent sources of HENs. In the more standard Waxman-Bahcall 
model, however, GeV emission, associated with higher Lorentz factor, results 
in decreased HEN emission (e.g. [30]). See e.g. [41] for a summary of possible 
emission mechanisms. 

Due to the small number of detected LAT GRBs, as well as our limited 
understanding of the emission mechanisms, investigating the temporal behav- 
ior of high-energy gamma-ray emission is difficult. It seems, however, that 
in many cases the light curves follow a power-law decay in time: Flat oc 
with a ^ 1.5 [51, 59, 60], while the rising parts for some cases are consis- 
tent with Flat oc t^. We note that the characteristic rising and decay were 
observed both for short and long bursts [51]. 

For the purposes of GW+HEN analysis, we are interested in the propor- 
tion of high-energy photons emitted in a given time window, T, following the 
onset of the GRB. For this estimate, following [52], we model the temporal 
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distribution of high-energy photons with a smoothly broken power law 



(4) 



with a = 2 and /3 = 1.5 as discussed above. The peak emission is at time 
tpeak = tb{a/ j3y^^°''^^\ A is an amplitude and B is the background level, 
both irrelevant for the calculation of the temporal behavior (for simplicity, 
we will assume B = 0). 

Due to the low number of available data, we restrict our analysis to esti- 
mating the ratio of GeV photons arriving in the 150s time window, i.e. the 
GRB emission time window described above in Section 3.1. Following [51], 
we use 8 out of the 12 GRBs for which the detected GeV emission has high 
enough signal-to-noise ratio for the reconstruction of the temporal structure. 
We assume that the emission of these GRBs follow Equation 4, with tpeak 
equal to the maximum shown in Figure 4. in [51]. Taking all 8 GRBs into 
account with equal weights, we verify that ~ 90% of the GeV photons arrive 
within a 150s time window after the onset of the GRB (see Figure 1). For 
comparison, 50% of the photons arrive within ~ 7s, while 95% of the photons 
arrive within 650s. We conclude that the 150s time window is a reasonable 
upper bound for the high-energy gamma-ray emission as well. We note that 
the obtained time window does not depend significantly on the initial rising 
part. The important factors are the peak time and the slope afterward. 

We note that, while the majority of GeV photons seems to arrive within 
the first 150s after the onset of the GRB, the importance of detected GeV 
photons from the perspective of neutrino emission can be that it might indi- 
cate a process, such as external shocks, that can produce HENs. If in such 
process the emission of GeV photons is not proportional to the emission of 
HENs (i.e. if the durations of the arrival of ~ 90% of GeV photons and HENs 
are different), the duration of HEN emission might be longer than 150s. 

3.3. Gravitational-wave Emission 

While gamma-ray and HEN emission from GRBs are practically indepen- 
dent from the specific mechanisms driving the (varying) relativistic outflow^, 
GW emission is closely connected to the central engine. 



^The distance from the central engine at which observable gamma-rays and HENs are 
produced is much larger than the central engine (e.g. [21]). 
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Various mechanisms can result in the emission of GW transients. Here we 
first consider stellar core collapses, the likely progenitors of at least some long 
GRBs. For these the collapse of the rotating iron core, as well as the follow- 
ing rebound in the inner core instantly produce a GW burst lasting for tens 
of milliseconds (see e.g. [25]). Other mechanisms following core bounce can 
also emit GWs shortly, within a few seconds after the collapse. These pro- 
cesses include rotational non-axisymmetric instabilities, post-bounce convec- 
tion, standing-accretion-shock instability (SASI), or non-radial proto-neutron 
star pulsation [25]. 

The collapse results in the formation of a protoneutron star, that can 
further collapse to form a stellar-mass black hole (e.g. [61] and references 
therein). The rotating, collapsing core might also fragment into two or more 
compact objects, whose coalescence can emit strong GWs [62]. After the 
formation of a black hole, weaker GW emission is possible from accreting 
matter around this central black hole that can fragment due to gravitational 
instability. The resulting black hole - fragment binaries can emit GWs de- 
tectable by Earth-based GW interferometers [63, 64, 65]. For high enough 
accretion disk masses, the self-gravity of the disk can result in gravitational 
instabilities such as spiral arm or bar formation, emitting GWs (e.g. [66], 
[67] and references therein). Non-uniform, non-axisymmetric accretion can 
also deform the central black hole that will emit gravitational waves during 
its ring-down [67]. Such GW emission processes can last for the duration of 
the accretion, i.e. possibly the duration of the GRB. 

Compact binary coalescence, the likely progenitor of short GRBs, is an- 
ticipated to be strong GW emitters in the sensitive frequency band of Earth- 
based GW detectors (e.g [37, 68, 67]). Double neutron stars, black hole- 
neutron star binaries and black hole-white dwarf binaries are considered as 
possible sources (e.g. [67]). Most of the GW output of such sources is emitted 
in the form of a very short, ~ 1ms long transient (e.g. [69, 67]). Furthermore, 
many short GRBs occur at low redshifts [70], increasing their significance as 
potential GW (and HEN) sources. 

Based on the models discussed above, GRBs are expected to be strong 
emitters of GWs. While core-collapse and binary merger are expected to emit 
the strongest GW transient, emission is possible during the entire duration 
when the central engine is active (see Figure 1). GW emission can consist 
of either one or more short burst (e.g. due to core collapse or the infalls of 
fragments of the accretion disk [63, 64, 65]), or longer duration GWs (e.g. 
due to the emission of rotational energy [66]). If the central engine consists 
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of an accretion disk, the activity of the central engine coincides with the 
emission of GWs. 

3.4- High-energy Neutrino Production 

In a GRB, HENs and gamma-rays are expected to be produced due 
to the variabihty of the central engine's activity, that results in fluctua- 
tions of the relativistic outflow, creating internal shocks in the ejecta (e.g. 
[71, 72, 32, 73]). These internal shocks accelerate electrons and protons in the 
outflow through the process of Fermi acceleration. Shock-accelerated elec- 
trons radiate their energy through synchrotron or inverse-Compton radiation 
(e.g. [72, 74]), emitting gamma-rays. Shock-accelerated protons interact with 
gamma-rays [pj) as well as with other, non-relativistic protons (pp), produc- 
ing charged pions and kaons [75]. Pions and kaons from these processes decay 
into HENs through [76] 

7r±,ir±^/i± + z/^(z7^) (5) 

The resulting muons from these interactions are expected to immediately 
undergo radiative cooling, making this secondary HEN production insignifi- 
cant [76, 77, 27]. For the proton energy range Ep ^ 10^ — 10^'^ GeV the pj 
process dominates HEN production, while outside of this energy range the 
pp process is dominant [75]. 

Since internal shocks in the relativistic outflow result in both gamma-ray 
and HEN emission, HENs are expected to be produced during the emission 
of gamma-rays (see Figure 1). 

For the case of long GRBs created by core-collapse supernovae, internal 
shocks in the relativistic outflow can occur even before the outflow emerged 
from the stellar envelope, therefore HENs are also expected to be produced 
before observable gamma-ray emission [27, 78]. While gamma rays emitted 
in this early phase cannot escape from the star due to its large optical depth, 
neutrinos have much longer mean free paths and may therefore pass through 
the stellar envelope. This pre-GRB neutrino emission is expected to precede 
the start of gamma-ray emission by up to 100s [27, 78] (see Figure 1). 

GRB afterglows are not considered here as a part of the coincidence time 
window. Afterglows are produced by the relativistic jet driven into the sur- 
rounding medium (e.g. [74]). The observed radiation, similarly to the case for 
prompt gamma-ray emission, is produced by synchrotron emission of shock- 
accelerated electrons. The energy distribution of protons is expected to be 
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similar to that of electrons [74], therefore the softer emission spectrum of 
afterglows indicate that protons might also be of lower energy. These pro- 
tons would need higher energy gamma-rays to produce neutrinos^, that are 
scarcely present in the afterglow spectrum. 

We note that, while HEN production in GRB afterglows is likely not 
significant, a few ultra high-energy neutrinos (UHENs) of energies ~ 10® — 
lO^GeV might be emitted during the afterglow phase [79, 74, 80], if GRBs 
can accelerate some protons to energies ~ lO^^GeV. In our time window 
estimation, we do not take into account such emission from the afterglow. 

3.5. GRB Precursors 

Gamma-ray bursts are sometimes preceded by fainter, softer electromag- 
netic emissions, so-called precursors (e.g. [81]). The underlying mecha- 
nism(s) creating these precursors is yet unknown. They have been detected 
for about 8 - 20% of GRBs [81, 82, 46, 83], with some GRBs apparently 
having multiple precursors. It is expected, however, that a higher percentage 
of GRBs are preceded by precursors, as many of these might be missed due 
to, e.g., beaming, low signal-to-noise ratio, the proximity of the precursor 
and the main event, or the definition of what is considered a precursor [81]. 

Analyzing precursors detected by BATSE, Lazzati [81] finds that the 
net count of a precursor, a value connected to the energy of the events, 
is up to about 1% of the main burst, most precursors having non-thermal 
spectra. Burlon et ai, using a more detailed analysis of data from the SWIFT 
telescope [82] and BATSE [46], finds that precursors, on average, emit ~ 
30% (for SWIFT) and ~ 10 - 20% (for BATSE) of the energy of the main 
burst. Burlon et al. find that precursors and main events have very similar 
spectral properties, concluding that precursors and main GRBs are hkely to 
be produced by the same underlying mechanism. 

There are several different GRB precursor models in the literature (e.g. 
[46]). The "two-step engine" model [84] argues that precursors are produced 
by initial weak jets, that are created by and coincide with the actual core col- 
lapse. The "progenitor precursor" model (e.g., [85]) connects precursors with 



•^In photomeson interactions producing HENs, the photon's energy and the proton's 
energy Cp are related, at the threshold of the A-resonance, by 

e^ep = 0.2GeV^r2 (6) 

in the observer frame, where V is the bulk Lorentz factor of the outflow [74]. 
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isotropic emission marking the jet breakout. The "fireball precursor" model 
[86] considers a simple model of an isotropic, post-acceleration GRB with 
constant rest mass and kinetic energy. Depending on the considered model, 
either or both the precursor and main GRB can have GW and/or HEN emis- 
sion. We note that, while the models above are concerned with precursors 
from stellar core collapses, precursor emission has also been identified from 
short GRBs [83]. 

We estimate an upper bound for the time difference between the onset 
of a precursor and the onset of the main burst using the results of Burlon 
et al. [46]. Burlon et al. analyzed 2121 BATSE GRBs, out of which 264 
(12.5%) was found to be preceded by one or more precursors (in total 369 
precursors). From this data, we determine an upper bound ^p^^^"'^''"'" to the 
time difference between the start of the precursor and the main burst, which 
is the 95% quantile of precursor time differences measured between the onset 
of the precursor and of the main GRB, as calculated by Burlon et al. Figure 
2 shows the distribution of the time differences and the 95% quantile, for 
which we obtain (also see Figure 1) 

^precursor ^ ^ 

For comparison, we note here that fP^^^'^'^^'^°^ ~ X50s an^l fP^'^'^^^'^"^ = t^^^^"^^"^ ~ 
350s. We use fP^^^'^'^^^°^ the upper limit of time difference between the main 
burst and a precursor. 

For the main GRB above, we used an upper limit of ~ 100s for the time 
difference between the activation of the central engine and the onset of ob- 
servable gamma-ray emission. We will use the same delay upper limit for 
precursors, i.e. the delay of the onset of gamma/X-ray emission from the 
precursor is estimated to be less than 100s after the start of the central en- 
gine (see Figure 1). 
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